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Abstract The large-scale structure of the magnetic field in the solar corona 
provides the energy to power large-scale solar eruptive events. Our physical 
understanding of this structure, and hence our ability to predict these events, is 
limited by the type of data currently available. It is shown that the multifractal 
spectrum is a powerful tool to study this structure, by providing a physical 
connection between the details of photospheric magnetic gradients and current 
density at all size scales. This uses concepts associated with geometric measure 
theory and the theory of weakly differentiable functions to compare Ampere’s 
law to the wavelet-transform modulus maximum method. The Holder exponent 
provides a direct measure of the rate of change of current density across spatial 
size scales. As this measure is independent of many features of the data (pixel 
resolution, data size, data type, presence of quiet-Sun data), it provides a unique 
approach to studying magnetic-field complexity and hence a potentially powerful 
tool for a statistical prediction of solar-flare activity. Three specific predictions 
are provided to test this theory: the multifractal spectra will not be dependent 
on the data type or quality; quiet-Sun gradients will not persist with time; 
structures with large current densities at large size scale will be the source of 
energy storage for solar eruptive events. 
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1. Introduction 

Solar flares and coronal mass ejections occur in active regions in the solar corona; 
however these magnetic structures are anchored in the solar photosphere. As 
magnetic flux emerges and is shuffled around the photosphere, the magnetic 
structure of the coronal field is stressed into non-potential configurations. Energy 
released in solar eruptive events must be present in the coronal field prior to the 
onset of the solar flare and/or coronal mass ejection, probably built up over hours 
to days as the active region either grows (by flux emergence) or as footpoints are 
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shuffled to result in sheared configurations. One of the key challenges in solar 
information processing is to use data of the photospheric magnetic field to infer 
the non-potential magnetic energy inherent in the coronal magnetic structure. 
The precise coronal conditions required to create these energetic events are as 
yet unknown, as coronal magnetography remains in its infancy. As we lack the 
ultimate dataset of 3D velocity and magnetic-field information throughout an 
active region, the prediction of solar flares remains a statistical problem 

One approach to this statistical problem is to use the photospheric magnetic- 
field measurement as the boundary for a 3D extrapolation of coronal structure 
( DeRosa et ah . 2009h . This yields important insights into our knowledge of 
coronal magnetic fields, one of which is that the final extrapolation is very 
sensitive to the boundary conditions {i.e. the spatial resolution, noise levels, 
image size, and type of magnetic-field data). As the field structure required to 
store enough energy to power a solar flare must be non-potential, there is often no 
unique coronal solution. Further the extrapolation is computationally expensive, 
the required vector data (even in the photosphere) are historically rare, and the 
chromospheric-field data (where the field may be force-free) remains tantalizingly 
out of reach. The best datasets in terms of sample size and homogeneity are full 
disk {i.e. low spatial resolution) and low-cadence longitudinal field data. The 
main source of these data are the Global Oscillations Network Group (GONG) 
and the Michelson Doppler Imager (MDI) onboard the Solar and Heliospherie 
Observatory (SoHO). Because of this data availability, a lot of research has 
focused on searching these data for a statistical measure of the coronal free 
energy from photosphe ric proxies ( Schriiver . 2007 : Leka and Barnesl . 2007 . 2008 : 


iBloomfield et al. , 2012 ) by correlating these proxies with solar flares as observed 


by GOES. This approach has displayed some potential in terms of a statistical 
solar-flare prediction, but often these proxies are not well physically motivated 
and are not well connected to theory. As such, it is difficult to realize the precise 
conditions under which these proxies fail {i.e. where they predict a flare that 
does not occur, a nd where they fail t o predict a flare) and what size scale is 
most appropriate ( Ireland et al. . 2008h . 


The multifractal spectra may provide a useful proxy of the build-up of mag¬ 
netic energy for in solar active regions. The proposition often advanced for this 
connection is that the scaling of magnetic structures is associated with the 
regime of fully developed turbulence at high Reynolds number. It is also based 
on the close analytic al connection between fractal geometry a nd turbulence i n 
magnetized plasmas ( McAteer. Gallagher and Conlonl . l2010l : McAteeil. 201^. 
Emp irical studies have used fractals on a large MDI dataset (iMcAteer et al 


2005 ). multifractals on a smaller MDI dataset (Conlon et al.. 20081. l2010l) . the 


multiscale Eourier index on a small MDI dataset (|Howet1^_et ajJ^ 20081 ). the 
standard Fourier scaling index on a large dataset ( Abramenkol. l2005l). and the 
struc ture function approach on a small data s et (Abramenko and Yiirchvshvnl . 
2010l ). However recent studies bv iGeorgoulisI ( 20I3I) and lMcAteer et al. ( 20151 ) 


reach a contrary conclusion that studies in multifractals do not provide a predic¬ 
tive ability for the onset of solar flares. One potential reason for this discrepancy 
is that multifractal scaling can also be produced by simple statistical similarity 
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in data of insufficient sp atial resolution, and so is not necessarily an indication 


of increased complexity ( Makarenko et al. , 2012ll . 


This article seeks to close the gap between these contradictory findings by pro¬ 
viding a complete analytical connection between multifractal formalism and the 
set of 3D equations that govern the small-scale and large-scale magnetic structure 
on the Sun. In Section [2] the background on the appropriate use of multifractals 
is provided with respect to the frozen-in field structure, arising from the large 
magnetic Reynolds number in the solar photosphere. Section [3] introduces the 
concept of multifractals and discusses the wavelet-transform modulus maxima 
approach of determining this spectrum from 2D images. The multifractal spectra 
is then connected to the magnetic structure as described by Ampere’s law in 
Section 21 Finally Section [5] summarizes this approach, discusses why previous 
results are contradictory, and makes specific predictions of this theory that can 
be tested empirically in order to make future advances in our understanding of 
how magnetic energy is stored in the solar corona. 


2. The Magnetic Reynolds Number 


It seems that the magnetic-field structure in active regions is the only viable 
means of storing and releasing sufficient energy to power solar eruptive events. 
A simple discussion of the equations used to describe this magnetic field leads to 
a couple of possible choices of characterizing the field. From first principles, the 
induction equation is given by eliminating the electric field between Faraday’s 
law and Ohm’s law 

^ = Vx{vxB)-W^, ( 1 ) 

(7 

which further reduces by Ampere’s law to 

f) Fi 

— = y X {vx B)-y X {r]S/ X B), (2) 

where rj = is the magnetic diffusivity. By Gauss’ law, and by means of 

a simple vector identity (and assuming constant ry), this further reduces to the 
familiar form of the induction equation, 

f) B 

— = V X {vx B)+r]y^B. (3) 


This states that any local change in the magnetic field is due to a combination of a 
convection term and a diffusion term. The battle between these two terms defines 
whether an active region will grow or shrink with time. The diffusion occurs 
by buffeting from supergranular flows, and the advection occurs by smaller- 
scale shuffling (granulation) or flux emergence. The ratio of these two terms is 
described by the magnetic Reynolds number, 


V X (v X R) 
r]V^B 


vl 

5 

V 


( 4 ) 
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which acts as an indication of the coupling between the plasma flow field and the 
magnetic field. For typical photospheric values in solar active regions (/ ~ 10® 
Mm, V « 10 ms“®, 77 ss 10® m^ s ®), the magnetic Reynolds number is much 
greater than one. In this ^ 1 regime, magnetic-field flux is advected with the 
plasma flow, until such time that gradients are concentrated into short enough 
length scales that diffusion can balance convection {I Ri 100 m). Ohmic dissipa¬ 
tion becomes important, and magnetic reconnection can occur. Essentially the 
large Rm allows for the build-up of energy, followed by a sudden release. Physical 
system with large i?ni are fully turbulent magnetic systems, and this naturally 
leads to studies of turbulence and complexity as potential key parameters to 
developing a deeper understanding of the physics behind active region evolution 
and energy storage required for solar eruptive events. 


3. The Multi fractal Spectrum 


The (mono-) fractal dimension of any object is a measure of the self-similarity 
across all size scales, or the scaling index of any length to area measure, [A oc P], 
where a is the capacity dimension, often just called the “fractal” dimension. In 
the more general case, a multifractal system will contain a spectrum of frac¬ 
tal indices of diffe rent powers f f(a)] (although strictly this is only true in the 
asymptotic limit, ( Turiel et al.l . boosf p. and takes account of the fractal measure 
at each point in space. The /(a) curve is equivalent to the D{h) multifractal 
spectrum, where D is the Haussdorff dimension (or filling dimension) of each 
Holder exponent in the data. If the system only has one Holder exponent {i.e. is 
a true monofractal) then this D{h) “spectrum” will be just one point. If a system 
contains a large number of sets of structures with dif fering Holder e xpon ents (i.e. 
is a multifractal), then this will be a concave curve. IConlon et ah ( 201 (ll l showed 
that flare-productive active regions have a multifract al spectrum peaking at large 
D {D > 1.2), agre eing with the monofrac tal study of McAteer et al. ( 2005ll . and 
high h (h > 0.3). Kestener et al.l ( 201fll l showed that it is possible to use this 
to separate quiet-Sun gradients from active-region gradients using the wavelet - 


transform modulus maxima (WTMM). The WTMM method (IMuzv et al.l . 119911 ) 


provides a novel means to calculate the multifractal spectrum of any dataset. 
A full description of the method as applied in ID is available in McAteer et 
al. ( 2007l l and iMcAteer and Bloomfield~(2013ll, as a pplied to solar flare hard 
X-ray data, and as applied in 2D by Kestener et al. ( 2010l l and McAteer et al 


( 201fll l. Briefly, this is achieved by initially computing the 2D continuous wavelet 
transform [T] of an image [f{x,y)]\ 

T^[/](b,a) = V(r,[/](b,a))=V(V*/) , (5) 

at each location [b] and where "0 is a smoothing kernel given separably into 


01 = 00/dx 
02 = OBjOy. 


( 6 ) 


This amounts to smoothing the data by a Gaussian kernel and applying a gradi¬ 
ent filter to the smoothed image to obtain T. The image is then smoothed by a 
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larger Gaussian kernel and the operation is repeated. Figure [T] shows gradients 
at three scales overplotted on magnetic-field data. Note that both quiet-Sun and 
active-region gradients are detected, as are the “fake” gradients at the center 
of the large negative spot (where the field incorrectly goes to 0 in the center 
of the largest spot). The largest gradients are at the mixed polarity “delta” 
components at the smallest scale. At the larger scales, the gradients tend to 
blend together into larger structures over the large-scale polarity-inversion lines 
evident in the active region. At the largest scale, edge effects creep in as vertical 
and horizontal streaks. At each scale, edges are extracted as maxima chains 
where the wavelet-transform modulus [M^[/](b,a) = |T.^[/](b, a)|] is locally 
maximum in the direction of T^[/](b, a). At each scale, the local maximum 
is then detected along each chain. A set of maxima lines, [Axo] is formed by 
connecting these points across scale. Statistically, these lines contain all of the 
information about the local Holder exponents [h] in the image. Along any one 
maxima line that connects down to any maxima xO in the limit at the smallest 
scale, the wavelet-transform modulus behaves as a power law with exponent 
h(xO): 

MMLM] (7) 


By tracking each gradient up through scale, from each location [xO] (where some 
gradient occurs at the smallest size scale), each of these locations can be char¬ 
acterized by its own Holder exponent that prescribes how that gradient varies 
across scale. The data can then be characterized as a histogram of these /i-values. 
A partition-function techniq ue is then used to con vert this histogram to the full 
D{h) multifractal spectrum ( Kestener et al. . 201fll l and provide appropriate error 
bars. 


4. Connecting the Multifractal Spectrum to Magnetic Structure 


The source of magnetic structure in the solar atmosphere is provided by Ampere’s^ 
law, 

V X = /tq J. (8) 

Where the constant is neglected, this breaks down to a combination of spatial 
gradients, 


V X B = 


i j k 
d/dx d/dy d/dz 

Bx By Bz 


{dB,/dy - dBy/dz) i+ 
{dB^jdz - dB^/dx)j+ 
{dBy/dx - dB^jdy) k 


(9) 


Where we are limited to measurements of B^ {i.e. longitudinal-field measure¬ 
ments near disk center), we neglect the other spatial gradients involving B^ and 
By. This is a good approximation near active regions with strong gradients. 
Assuming flux [/ B.da] is conserved, then as density falls off exponentially with 
height, a flux tube will expand in area and so decrease in field strength. Hence 
d/dz will fall off rapidly, with the same exponential form as density. Neglecting 


SOLA: ms.tex; 29 June 2015; 0:31; p. 5 










R.T.J. McAteer 


spatial gradients involving and By may not turn out to be a viable assumption 
in some cases (something to be accounted for in handling full magnetic vector 
data). Proceeding with B^ only, Equations ([5]) and (O become 

J{x,y) = {dB^/dy)i-{dB^/dx)] , (10) 

and hence 

\J\ = ^{dBjdyf + (dB^/dxf = V^^yB, . (11) 

Calculating magnetic spatial gradients is a non-trivial issue. Gradients may occur 
in any direction, and are not limited to the pixel resolution of the data. The most 
appropriate way to measure gradients is to do so for every angle and every spatial 
size scale, and then to study how these vary over size scale, [s], i.e. 

d\J\/ds = dV,,yB,lds. (12) 

By adopting a derivative of Gaussian as the mother wavelet in the analysis of 
WTMM in Equations ([S]) and (jB]), the local Holder exponent, h(xO) in Equa¬ 
tion © is the rate of change of the gradient of the magnetic field with size scale, 
i.e. in the asymptotic limit. 


k(xO) = dVx,yB;^/ds. 


(13) 


For the gradients identihed in Figure [TJ the modulus maxima are identified 
along each ridge, and then tracked up through scale as described in Section [S] 
Figure [5] shows the result of this tracking procedure. Each gradient is tracked 
to progressively larger and larger spatial size scales. According to Equation ([T]) 
the slope of each of these lines is the Holder exponent of each gradient tracked. 
According to Equation (IT^ , this value is also the spatial rate of change of current 
density. 

The lines in Figure [2] are color-coded according to the h value. Glearly some 
gradients exhibi t a negative or nea r-zero h. Positive, near-zero Holder exponents 
are “dust-like” (iGeorgoulisl. l2002h , meaning they are sparse and only exist at 
small scales {i.e. the power drops off quickly with increasing scale). By defi- 
nition in Equation (I71l. the Holder exponent is positive, however recent work 


(ILeonardiizzi et al.l . 120141 1 has extended this analysis into a negative regime via a 


secondary component. In Equation (1121) . both negative and near-zero exponents 
correspond to locations where some current density may be present at small 
scales, but not at large scales. It is clear from Figure [3] that these locations are 
in the quiet Sun. When tracked over time, these features will not persist and 
hence they cannot build up the magnetic energy required to be the source of 
solar flares. 

The positive Holder exponents (blue) are large positive h. These are “strong”, 
meaning they exist at small scales and are larger at larger scales. In Equa¬ 
tion (nm, these are locations where current density increases at large scale. It 
is clear from Figure [3] that exponents are preferentially located in and around 
strong magnetic fields. In particular, there are groups of strong h (blue) to the 
North and South of the active region between the strong negative spot and the 


SOLA: ms.tex; 29 June 2015; 0:31; p. 6 















Frozen-in Fractals 


positive regions. There is another group of strong h (blue) to the West (right) at 
the negative polarity plage. All the real active-region gradients present an Holder 
exponent greater than 0.7. The expectation is that when tracked over time, these 
are the gradients that persist to provide magnetic structure for strong current 
densities that can exist up to larger spatial scales. It is expected that these are 
the locations where magnetic energy can therefore build up over time and these 
will be the locations of energy release. 

Figure H] (left) shows the value of each Holder exponent plotted against the 
gradient at the smallest size scale. Interestingly, a large gradient at small size 
scale is not necessarily an indicator of a large Holder exponent. This shows 
the need to consider gradients across scale, rather than just gradients on a 
pixel-by-pixel basis. The frequency histogram of these Holder exponents (over 
plotted) is distinctly bimodal. There is a small negative-/i distribution and a 
larger, long-tailed positive-/i distribution. Figure S] (right) shows the same plot 
of h-values, but this time overplotted with the multifractal spectrum preselected 
only for positive /i-values that show a large gradient at small size scale, i.e. 
this multifractal spectrum is formed from only those data in Figure [2] that are 
blue and show a gradient larger than 200G px“^ at the smallest size scale. The 
multifractal spectrum is calculated using a partition func tion algorithm a pplied 
to t he gradients of Figure [2] and described in detail bv iMcAteer et ah J2007 1 
and iKestener et iah ( 20ldl . This technique of preselecting Holder exponents is 
justified as the quiet Sun presents its own individual multifractal spectrum; 
therefore if all the Holder exponents are used, the resulting (averaged) multi¬ 
fractal spectrum would be incorrectly broad and would incorrectly peak around 
h = 0. By preselecting on both gradient size and Holder exponent it is possible 
to by remove quiet-Sun contamination and therefore extract the active-region 
multifractal spectra from any dataset. 


5. Conclusions and Future Work 


The multifractal spectrum exhibits much potential in providing a suitable, and 
physically mot ivated, proxy f o r the complexity (and hence non-potentiality) of 
active regions. Ireland et al.l ( 2008h showed that a multiscale neutral-line de¬ 
tector could be useful for studying solar flares. Together with the fact that the 
magnetic Reynold’s number is very large, the application of the multifractal 
spectrum is appropriate. As Ampere’s law breaks down into a collection of 
magnetic spatial gradients, the route to obtaining the multifractal spectrum 
through using these gradients is optimum. By comparing Equation (1101) for the 
horizontal components of the current density with Equation o that defines 
the Holder exponents it is shown that the Holder exponent explicitly defines 
the rate of change of current density with spatial size. The histogram of these 
Holder exponents then explicitly defines the overall magnetic complexity of the 
image. It is important to realize that the histogram in Figure |4] is distinctly 
bimod al. This bimodality also agrees wi t h pre vious work on solar magnetic 
helds ( Kestener et al.l . [20101 : IConlon et al.1 12 OIOI I and this bimodality is hard to 
capture in the usual Legendre-transform approach to the multifractal spectrum. 
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It only becomes evident when each gradient is tracked and studied separately. 
It can represent a phase transition between two distinct measures. Clearly this 
bimodality must be studied further across a la rger sample of mag netic-field data. 
This bimodality e xplains why re c ent research ( Georgouli^. 2013 ) fa iled to repeat 
earlier findings of IConlon et al.l ( 20081 . 2010l l and Kestener et al. ( 2010l 'l. If the 
entire histogram is incorrectly fitted with a single distribution, the resulting 
incorrect multifractal spectrum will be wide and will peak at lower D and lower 
h. This is a direct result of the influence of the quiet-Sun gradients influencing 
the overall result. Further, as quiet-Sun gradients are the weakest, it is easy to 
lose these from image to image, or from instrument to instrument - hence a cross 
comparison or study of any time sequence of data that incorrectly includes these 
quiet-Sun gradients will result in inconsistent, fluctuating, multifractal spectra. 

It is important to note that the measurement that is taken from the data is 
that of how the gradient changes with size scale - this forms the first prediction 
that the same solution will be determined irrespective of the data type. Future 
work must test this prediction by comparing the multifractal spectrum on co¬ 
temporal and co-spatial regions across datasets from multiple instruments. It 
is this potential to provide a link between the large MDI dataset, the large 
GONG dataset, and all future datasets that makes this prediction so intriguing 
for statistical solar-flare prediction. A second prediction arising from the work 
presented here is that those gradients defined as quiet Sun - those red and 
green points in Figure [2] and Figure [3] - should not persist with time. These 
gradients can provide some direct heating to the solar atmosphere, but they 
cannot persist, and hence energy cannot build up to large values of current 
density. This prediction can be tested by tracking these from image to image 
- the gradients defined as quiet Sun will persist for much less time than those 
gradients defined as active region. A separate study of negative exponents in 
the quiet Sun ma y provide insight into these small-scale magnetic flux densities 
(see discussion by Gharon and Trou^ ( 201,‘lli '). The third prediction is that the 
energy required to power solar eruptive events comes from those large gradients 
with large Holder exponents. Larger solar flares should result from the largest 
gradients that persist for the longest time - indeed it may be possible to calibrate 
the actual free energy stored in the magnetic field by comparing the size and 
persistence of both the overall multifractal spectra and the individual gradients 
with the resulting energy releases. If it can be shown that the devil of the largest 
scale structures is provided by the detail of the smallest scale gradients, it will 
then be possible to combine large-scale magnetic proxy studies with smaller-scale 
extrapolation work to provide future insight into coronal structure, a study that 
can then be further tested when coronal magnetic-field data become readily 
available. 
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Figure 1. A multiscale gradient decomposition of an active region magnetogram. The back¬ 
ground image in each row is the same magnetogram, covering an area of 540Mm by 270Mm, 
scaled to -h/- 1000 G of the vertical color bar. Top: Contours of magnetic gradient (G/Mm) 
according to the horizontal colorbar, and calculated a spatial scale of 0.6 Mm. Middle: Contours 
of magnetic gradient on same color scale as top, but calculated at a spatial scale of 1.7 Mm. 
Bottom: Contours of magnetic gradient on same color scale as top, but calculated at a spatial 
scale of 6.3 Mm. 
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Figure 2. The Holder exponents of the data presented in Figure^ Each gradient is tracked to 
larger and larger scale. The value of the gradient at each size scale is plotted in log-log space. 
The resulting slope of each line is the Holder exponent of each gradient and describes the rate 
of (spatial) change of current density fEauation I12I I. Red lines show a negative value for h, 
green lines shows values of h around zero, and the blue lines show positive h values where the 
current density increase for larger scales. The same vertical colorbar is used in Figure 0 
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Figure 3. The location of each gradient is described by its Holder exponent. The negative 
(red) and small (green) values of gradient are mostly in regions of quiet Sun. The strong 
positive (blue) Holder exponents are the locations where strong gradients exist up to large size 
scales. These blue gradients delineate the polarity inversion lines within the active region. 
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Figure 4. (Left: A scatter plot of Holder exponent against the gradient at smallest size scale. 
The histogram of this entire distribution peaks at around h = 0. Right: The multifractal 
spectrum is calculated only for those points that have small-scale gradients greater than 
200 G/px, horizontal line plotted on the left, and with positive Holder exponent. The resulting 
multifractal spectrum peaks around D{h = 0.8) = 1.35. Error bars are only plotted for those 
values where the calculated error is larger than the diamond symbol. 
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